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Идентификация параметров математических 
моделей динамических систем управления 


В статье рассматривается программный комплекс для оценивания параметров математических моделей 
сложных динамических систем на основе экспериментальных данных, полученных в процессе их 
функционирования. Математической базой данной системы являются современные методы оптимизации, 
рекуррентные алгоритмы оценивания, методы линейной алгебры, методы линейного и нелинейного 
регрессионного анализа, методы вычислительной математики. Комплекс использует объектно- 
ориентированный подход при реализации программ, имеет удобный интерфейс, не требует от пользователя 
знания особенностей языка программирования или деталей схемы организации процесса вычислений в 
комплексе. 


Введение 


На современном этапе развития производства и конструирования сложных 
динамических объектов невозможно обойтись без предварительной разработки ма- 
тематической модели такого объекта. Для тестовых испытаний строятся сложные 
математические модели (что, как правило, задаются системами дифференциальных 
уравнений). Такие модели на практике очень сложные и содержат много параметров 
и параметр управления, с помощью изменения которого возможно тестирование 
модели при различных условиях. Ясно, что важнейшая задача такого моделирования 
состоит в том, чтобы как можно точнее описать модель. Но не всегда при описании 
моделей можно определить значение некоторых параметров. Примером таких величин 
могут быть внешние силы природы, действующие на поверхность самолета, внутренние 
волновые потоки. То есть те величины, которые не может контролировать человек. 
Одной из реальных задач идентификации была идентификация параметров модели 
продольного движения самолета АН [1]. Практически можно было исследовать влияние 
рулевого крыла на динамику полета самолета и моделировать движение и при этом 
находить такие параметры, как аэродинамический момент, подъемную силу, 
сопротивление внешней среды. 

Следовательно, задача заключается в том, чтобы по экспериментальным данным 
приблизить математическую модель к реальной с помощью критерия. Как показывает 
практика, методы минимизации такого критерия зависят от структуры самого критерия, 
от общего вида системы. Итак, для моделирования и исследований необходимо иметь 
достаточно широкий спектр методов минимизации, имеющих различные характери- 
стики скорости сходимости. Каждый из существующих алгоритмов многомерной 
оптимизации имеет свои преимущества и недостатки. Поэтому при решении задачи 
минимизации функции был использован разработанный комплекс алгоритмов мини- 
мизации, последовательно подключающий различные алгоритмы на тех этапах процесса 
оптимизации, где каждый из них является наиболее эффективным. Существенной помо- 
щью в выборе эффективного метода было графическое изображение функции-критерия. 

Как видно, задача идентификации приведена к задаче минимизации, что является 
само по себе достаточно сложной задачей. Но на данном этапе возникает проблема — 
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достаточно неудобно было бы для каждой новой системы или нового метода мини- 
мизации перекомпилировать проект. Цель работы заключается в том, чтобы создать 
такой комплекс алгоритмов, который бы легко адаптировался под разные системы и 
разные методы. Каждый пользователь может написать свой метод или описать 
новую систему, оформить в виде О, и тогда автоматически этот метод или система 
будет заметна программе. Естественно, все методы и системы должны подчиняться 
соответствующему шаблону. Таким шаблоном выступают базовые классы, от которых 
пользователь должен наследовать свой класс. Для примеров и тестирования некоторых 
систем реализованы методы одномерной, многомерной минимизации. 

Поскольку на практике мы имеем достаточно сложные системы, невозможно 
заранее прогнозировать, насколько корректным будет решение и какие методы необхо- 
димо использовать. Как видно, решением будут минимизированы параметры для 
функции критерия, поэтому понятно, что для того чтобы подобрать удачные методы 
минимизации нужно знать топологическую структуру функции критерия (ее графи- 
ческое представление). Поэтому была реализована возможность посмотреть на функцию 
в трехмерном пространстве, фиксируя другие переменные. Такая предварительная ин- 
формация позволяет увидеть, какая форма критерия — выпуклая, овражная или может 
существовать несколько экстремальных точек — по этой информации можно подобрать 
метод, как лучше решать данную задачу. 

С помощью данной системы был решен ряд тестовых примеров, которые показали 
ее работоспособность и удобство подключения новых систем. Также была взята 
реальная система, описывающая продольное движение летательного аппарата, и подсчи- 
таны параметры на реальных данных. 


Постановка задачи и критерии оценивания 


Рассмотрим постановку задачи идентификации для динамических систем, описы- 
ваемых обыкновенными дифференциальными уравнениями, правые части которых 
нелинейные относительно параметров и фазовых координат, а также проведем класси- 
фикацию задач и алгоритмов их решения в зависимости от классов модели и кри- 
териев идентификации. Запишем общий вид модели 


Ах 
ИО Р (1) 


Начальные условия для задачи Коши "(0 )=ху Вектор параметров 


рР=(р......Р»). Вектор модельных фазовых координат 


хе )= ( © . 


Вектор экспериментальных фазовых координат 
е е е 
2 = (40,0) 


На основе экспериментальных данных (наблюдений) за объектом на интервале 
времени [4 +(А-ПТ, & +АТ], К=1...,г на каждом из таких промежутков фикси- 


руются экспериментальные значения фазовых переменных, обозначим их как х” ‚ где 
а ВА Е.М МЕТИМ ум 


Необходимо идентифицировать (определить) параметры р из условия соответ- 
ствия выбранной модели (1) данному объекту по некоторому критерию качества модели 
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Гр) так, чтобы Г(р”’) = тт Г(р) ‚ при этом параметры р могут быть произвольными 
р 


или иметь какие-то ограничения реЕД={р:а, < р, <Ь}. 


В зависимости от требований, предъявляемых к задаче, выбирается соответству- 
ющий критерий качества /(р), который может задаваться разными способами: 


1. В виде функции среднеквадратичного критерия: 


и-+Т 
Кр)= | к” (р)- хр ‚ || | — есть расхождение модельными фазовыми ко- 
И 
ординатами и фазовыми координатами экспериментальных данных, || || - норма в 


Евклидовом пространстве, соответствующая ему дискретная форма 


№ п 
2 
Кр)= У У" е,,р)-*а,). 
а У У 
Д=Н= 
2. В виде минимаксной функции несогласованности экспериментальных и модель- 
НЫХ фазовых координат: 


11 Кр= шах |х”@.,р)-х’@.)|| — сумма по всем максимальным отклонениям; 
1<=7<=М№ 7 


п 
т | тах ху (1 р р-х (1 7 )] — максимум по всем максимальных от- 
1=1 1<=1<=М№ 
клонениям. Как видно, здесь использовано Чебышевскую норму. 
3. В виде функции-невязки, которая в случае линейной функции правых частей 
системы (1) относительно параметров Р позволяет свести задачу к линейной 


идентификации (в данной работе линейный случай рассматриваться не будет) 
2 
ах № п 
Кр) = А -ЕОеьр) , Кр)= У Уре )- Боба риьвр. 
4 о 1 о 
Л=НЕ 
Выбор конкретного критерия оценки /(р) зависит от характера данных, от класса 


модели и требований, предъявляемых к алгоритму. 
При выборе первого критерия фактически решается задача среднеквадратичной 


аппроксимации данных эксперимента х’(Р) решениями задачи Коши (1) и (2), то есть 
значениями х” (+, р). Для произвольных типов моделей (1) этот критерий является 


наиболее удобным с точки зрения различия х”(.р) и х’(®. Функция Г(р) этого 
критерия является выпуклой, что гарантирует хорошую сходимость некоторых 
методов минимизации. Эта функция дифференцирована, позволяет применять гра- 
диентные методы спуска, которые являются достаточно удобными в вычислении. Но при 
применении методов поиска минимума гарантировано большое количество вы- 
числений функции /(р), а в некоторых методах первой или второй производной, что 


делает часть этих методов совсем несостоятельными к решению таких задач, по- 
скольку вычисления этой функции для каждой компоненты вектора параметров р 


требует решения задач задачи Коши на интервале [1, & +7]. 


При выборе второго критерия фактически строится Чебышевское приближение 
(равномерное приближение) экспериментальных данных х’({) и результатов задачи 


Коши х”(Ьр), минимизирует максимальное отклонение модельных решений и 


экспериментальных наблюдений на интервале [1, 1 +Т]. В некоторых случаях 
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практическое использование этого критерия может быть не эффективным. Когда 
наблюдения х’ (Г) наблюдаются на фоне случайных возмущений с большими выбросами, 
в этом случае необходимо предварительное сглаживание экспериментальных данных. 
Методы решения задачи идентификации при этом критерии, а также трудности 
использования алгоритмов остаются те же. 

Второй проблемой этого критерия является проблема недифференцируемости 
функции /(р), что заставляет использовать различные методы к нахождению градиента 


функции. 
Если рассмотреть различные виды дискретной формы, то можно сказать, что 


п 
критерий типа /(р)= У [ шах ху" ((;,р)-х’(;)] является достаточно целесо- 
Л АЯ 
1=1 1<=1<=М№ 
образным, поскольку минимизируется сумма по всем максимальным отклонениям, а 
значит, гарантирует уменьшение максимальных отклонений для всех максимальных 
отклонений. 


Для функции /(р)= тах [ шах |х"(.,р)-х(.)]] берем максимум по всем 
1<1<п 1<}<№М 1 1 РЯ 

максимальным отклонениям, а значит, уменьшать будем только само максимальное 

отклонение. 

Третий критерий наиболее распространен в инженерной практике. При этом 
критерии фактически решается задача среднеквадратичной аппроксимации производных 
экспериментальных данных функциями правых частей. В случае линейной по пара- 
метрам модели этот критерий позволяет свести задачу нахождения параметров р к 
линейной задаче, то есть к системе линейных уравнений Ар =Ь, менее сложной за- 
даче. Для этого критерия также есть трудности, во-первых, задача дифференцирования 
экспериментальных данных, которые наблюдаются, — это некорректная задача, во-вторых, 
во многих случаях система Ар =Ь имеет плохо обусловленную матрицу, а значит, 
решения могут быть как угодно далекими от истинных. Использование специальных 
методов решения линейных систем, имеющих плохо обусловленную матрицу, требуют 
расчетов с более высокой точностью, тем самым подавляя точность эксперимен- 
тальных данных. 


Общая характеристика алгоритмов нелинейной 
идентификации 


Общий алгоритм решения поставленной задачи предусматривает такие этапы: 

1. Найти решения системы дифференциальных уравнений (1) для начального 
(априорного) значения параметра р=(р.,...,р,), для начальных условий задачи 
Коши (2) на промежутке [1,1 +7]. 

2. Для данных значений х”(,р) — решений задачи Коши, х° (Г) — экспери- 
ментальных данных найти значение функции — критерия качества /(р), которая 
была описана выше. 


3. Найти р = аго шт /(р), найденное р — будет следующим значением параметра 
РЕ" 


для решения задачи Коши. 
Поскольку задача состоит в нахождении неизвестных параметров нелинейной 
системы дифференциальных уравнений (1) как относительно самого параметра, так 
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и относительно фазовых переменных, то основной задачей является задача много- 
мерной минимизации. Сформулируем задачу безусловной минимизации функции многих 
переменных, которая выбрана из одного из критериев качества, приведенных выше. 

Для решения задачи минимизации нахождение точки минимума функции при- 
меняются методы, которые включают в себя методы дифференцирования (а вернее, 
нахождение градиента). Если описать вычислительную схему — большинство алго- 
ритмов можно выделить в ряд последовательных этапов, а именно: 


1) выбор начального приближения параметров Ро и вычисления /(ру); 


2) определение направления спуска Р;, для минимизации [(р) (способ решения 
этой подзадачи является определяющим для каждого конкретного алгоритма) характерна 
своя стратегия при выборе направления спуска р, фактически — это нахождение гра- 
диента функции /(р); 

3) определение точки минимума 1, на луче р=р,+1*й,*е,, где 
пит (РАЯ *Й, * 8, ) = КР + Я * №№ * 8+); 

4) переопределение полученной точки минимума на луче р,._р, +А*Й, *2, в 
качестве исходной на следующем этапе: р. =р,+1*й *в,. 


Условный переход на 2), если критерий выхода из алгоритма не выполняется. 
Однократное выполнение этапов 2 — 4 будем называть одной итерацией процесса 
минимизации. На каждом из этапов нужно фактически решить некоторую само- 
стоятельную подзадачу. Для решения подзадач 1 и 3 можно использовать несколько 
эквивалентных (в смысле решаемой нами задачи) методов, которые реализованы в 
виде самостоятельных программных модулей. 

Для отдельной задачи минимизации выбор начальной точки хо осуществляется 
управляющей программой минимизации. Если начальная точка не задана, то вместо 
нее в качестве исходной точки для поиска выбирается случайная точка, что 
равномерно распределена в некоторой заданной области Ру в п-мерном парал- 


лелепипеде РО =(х:а<х<Ь. 
Для случая критерия качества равномерного приближения мы имеем недиф- 
ференцированную функцию /(р) = „тах, тах | хг (#,,р)-х; (#,) |. Невозможно исполь- 


<=71<=М 1<=< 


зование прямых методов нахождения р = аго т /(р) с помощью градиентов, посколь- 
рек" 


ку функция имеет вид Ё(1) = тах | х/ (1, р)-х, (#,) | (рис. 1). 
1<=<=и ' / ' / р 


Е) 


Рисунок 1 — Функция Ё(И) = тах | х/(Е.,р)-х: (Е,) 
1<=1<=п | 1 р Е 


Здесь много экстремумов, следовательно, нахождение градиента (напра- 
вления спуска) для одного максимума не гарантирует спуск другого максимума. 
Поэтому решение проблемы заключается в том, чтобы выделить максимумы 


Е тах (1 ),..., Е тах п(#,) для каждого такого промежутка (обозначим, А-й проме- 
жуток), решить задачу Коши (1) на промежутке [1,1 +1]. Найти для этих решений 


(р), найти для этого промежутка свое направление спуска р, . 
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Нахождение направлений спуска (на каждом из таких интервалов) уже позволит 
применять градиентные методы, которые мы описали выше. Пусть нашли такие 


векторы направлений спуска для каждого из 5 интервалов: ру,...,р,. Найдем теперь 
для линейной оболочки, построенной на векторах: ру,...,р,, такой вектор х ‚ который 
будет принадлежать этой оболочке, и будет наименее удален от начала координат. 

Пусть для каждого промежутка нашли вектор наискорейшего спуска р; (для 1-го 
промежутка) строим вектор х, который является линейной комбинацией векторов р; 
х=Ар!+...-+К,Р,. 


Решаем задачу минимизации 
|х||= (х,х) = К? | Р1]] +2 * К1(рЬр2)+....+ А5' || р5 |-> ша, 


5 
если м =1. Для этой задачи строим функцию Лагранжа 
15 


(Я) = 2? | РИ +Е2 * А(РЬр2)+....+ 5? || ру || 


+1* (УК — 1? || РИ || +2 * А1(РЬ р2)+....+ 45° || 25|), 
1= 
ОГ(А,К1,..К5) 0. ОГ(А,К1,..А5) _ 
Ок, | 04 
По сути, это линейная система уравнений. Ее решения относительно % и дадут 
значение искомого вектора х. 


Итак, если есть векторы ру......, р» › то для построения вектора х =Ар! +...+А,р», 


0; 


: 
при условии, что |х| > шти Ук, =1. Если решить систему, получим 
т 


2 
й а. 
2 


и) 
о 
[р |р, 


Рассмотрим подробнее подзадачу 4. В ней понятие «критерий» включает в себя 
два условия: 

1) алгоритм «исчерпал себя», что отразилось на выполнении так называемого 
«внутреннего» критерия (этот критерий свой для каждого алгоритма), тогда выход из 
программы минимизации происходит с кодом возврата Ки =1. 

2) было выполнено заданное число К/Т итераций, МтМ - это экземпляр класса 
многомерной задачи минимизации, однако внутренний критерий не исполнился, тогда 
возврат происходит с Ки=2, число МММ -— >1ЙТ задается пользователем или по 
умолчанию выбирается управляющей программой. При возвращении в управляющую 
программу минимизации с кодом Ки’ = 2 осуществляется анализ скорости сходимости 
поискового алгоритма с помощью подпрограммы 5ап. 

Если сбор информации автоматизирован, то можно придерживаться таких 


оценок сходимости. Если не выполняются условия 


1/2 1/2 
МИ т МИ т 


ое | бы" 
р ГЕ 
В — 


= К=1 


х МИ ГА р МИ УК 
то скорость сходимости данного алгоритма неудовлетворительная и пользователь 
должен изменить методы минимизации. В случае неудовлетворительной СХОДИМОСТИ 
пользователь также может управлять алгоритмами, изменяя значения их входных 


параметров или конфигурацию модулей, решающих задачу. 
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Заключение 


Разработана усовершенствованная версия этой системы для персональных 
компьютеров в среде \т@4о\5$ на С++, в которой представлены более 15 методов 
оценивания параметров и состояний, позволяющих решать как линейные, так и нели- 
нейные задачи идентификации в условиях «зашумленности» и неполного поля изме- 
рений, включая возможность обработки данных в режиме реального времени. Возможна 
переориентация комплекса на решение задач идентификации для других предметных 
областей применения. Данная система позволяет решать задачи идентификации для ди- 
намических систем, математическая модель которых описывается системами 
обыкновенных дифференциальных уравнений или функциональными соотношени- 
ями. В системе предусмотрен удобный интерфейс с пользователем на всех этапах 
решения задачи, удобная графическая поддержка процесса решения задач. 
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Петрович В.М. 

дентифжац параметр математичних моделей динам!чних систем керування 

У стати! розглядаеться програмний комплекс для ошнювання параметрв математичних моделей складних 
динамичних систем на основ! експериментальних даних, отриманих в процес! 1х функцонування. 
Математичною базою дано! системи е сучасн! методи оптимзацй, рекурентн! алгоритми ошнювання, методи 
лийно! алгебри, методи лейного 1 нелеййного регресйного аналзу, методи обчислювально! математики. 
Комплекс використовуе об’ектно-орентований щдхд при реатзацй програм, мае зручний 1нтерфейс, не 
вимагае в1д користувача знання особливостей мови програмування або деталей схеми оргавзацй процесу 
обчислень в комплекс. 


Ретомсй И.М. 

Тдепийсайоп оЁ Рагатеег$ о Маетайса| Моде о? Оупапие Сопб`о1 Зу$вет$ 

Те агасе 415си5$е$ а зоЙ\таге расКазе ог езитайпе Фе рагатеегз оРтлаветлайса! плоде! оРсотар!ех дупаптса! 
зузетз Базе оп ехрегитепа дайа оМатед ш Фе соигзе оЁР Феи: орегайоп. Те плафеплайса! Ъаз1$ оЁ 1$ зу$ет 
аге фе тодегт тефо45 оР орниманоп, гесигяуе езштаноп а!еог плз, Ниеаг аеебга, плефо4$ оЁ Ппеаг апа 
попПпеаг гезтеззюп апа[уз1, пе о4$ оЁ сотприайопа| таетайс$. Сотрех и5ез обес-опете4 арргоасВ шт Фе 
пиретепеаноп оЁ ргоотапаз, Ваз а сопуешепе пщеасе фа 4оез поё гедине апу Кпо\ еде оЁ ргоогапиие 
1апопасе ог Фе де 1$ оРШе эсБетез оЁ сай1сайопз ш фе сотрех. 


Статья поступила в редакцию 06.06.2011. 
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